Trap Models and Slow Dynamics in Supercooled Liquids 
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The predictions of a class of phenomenological trap models of supercooled liquids are tested via 
computer simulation of a model glass-forming liquid. It is found that a model with a Gaussian 
distribution of trap energies provides a good description of the landscape dynamics, even at 
temperatures above T c , the critical temperature of mode-coupling theory. A scenario is discussed 
whereby deep traps are composed of collections of inherent structures above T c and single inherent 
structures below T c . Deviations from the simple Gaussian trap picture are quantified and discussed. 



Glass-forming systems are ubiquitous in nature, and 
constitute the basis for wide-ranging technological appli- 
cations [jl, 0j. Furthermore, the theoretical and computa- 
tional techniques developed for the study of glass-forming 
liquids have been extremely useful in the study of static 
and dynamic phenomena in other complex systems, such 
as spin- glasses, polymers and proteins [^]. Unfortunately, 
despite recent progress, the underlying microscopic origin 
of slow dynamics in supercooled liquids remains largely 
an unsolved problem. 

The mode-coupling theory (mct), pioneered by Gotze 
and co-workers, is unique in that quantitative predictions 
are made for the viscous slowing-down exhibited by dy- 
namical correlators with the input of structural informa- 
tion alone [Q. In this sense, mct is perhaps the only 
ab initio theory of slow dynamics in supercooled liquids. 
mct appears to be a quantitative theory of cage forma- 
tion at intermediate times, but has some shortcomings in 
describing cage relaxation at long times. In particular, 
mct overestimates the the cage effect, leading to a "glass- 
transition" at a temperature T c > T g where T g is the 
experimentally determined glass transition temperature. 
Theoretical attempts have been made to include "hop- 
ping" processes that restore ergodicity below T c , however 
such theories currently demand a level of approximation 
that reduces the ability to make quantitative predictions 
1- 

A seemingly different perspective is provided by the 
"landscape" picture of slow dynamics in supercooled liq- 
uids |^|, gj. Here, one attributes slow structural relaxation 
to the complex pathways that connect configurational 
states or "inherent structures" (is) on the multidimen- 
sional potential surface. Wales and coworkers have per- 
formed extensive studies of the details of the landscape 
properties of model glass forming liquids @, ||. In par- 
ticular, they have noted that the landscape of Lennard- 
Jones liquids has a multi-funnel structure. The inherent 
structures residing inside a funnel are separated by small 
barriers, while barriers separating funnels may be large. 
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Heuer and coworkers ||, [T(| |n) have recently shown that 
the phase point of an unbiased trajectory is trapped in a 
"meta-basin" for long periods of time, making frequent 
small hops within the meta-basin, and infrequent excur- 
sions from one meta-basin to another, a scenario that 
that resonates with the multi-funnel picture. 

On the other hand, recent studies of mean field p-spin 
glasses H have established a deep connection between the 
mct and landscape pictures. The mct transition tem- 
perature T c appears as the point below which the only 
available saddle points of the energy are minima, whereas 
above T c , the vast majority of saddles are higher order 
saddles p|. Within mct (or mean field p-spin glasses), 
the dynamics freezes at T c in the lowest available 'criti- 
cal' saddle that percolates in phase-space, but never pen- 
etrates the lower minima region below T c . This is related 
to the above mentioned inability of mct to describe hop- 
ping events. 

The relevance of this result for Lennard-Jones systems 
was discussed in the important work of Angelani et al. 
]l3t and Broderix, Cavagna et al. |l~4|, [l5|, where the 
critical mct temperature T c , as identified from the ex- 
trapolated divergence of dynamical time scales, coincides 
with the point where the number of unstable saddles van- 
ishes. However, due to the non mean-field nature of the 
system, several important differences with MCT appear. 
First, because barriers are finite, the dynamics below T c 
is now dominated by activated transitions between basins 
of local minima (is), and not by the exploration of the 
"critical" saddle. As was shown in |l7), there indeed 
exists a separation of time scales between vibrational mo- 
tion within an is, and hopping between is, as envisioned 
long ago by Goldstein ^j. Second, the vanishing of sad- 
dles below T c and the inexistence of minima above T c 
are no longer sharp statements. One of the primary goal 
of the present paper is to actually show that the long 
time dynamical properties are indeed governed by these 
rare, deep traps, even above T c . The landscape studies 
of Wales et al. J?], |§| and Heuer et al. ||, |lj suggest 
that an activated dynamics of a more complex variety, 
namely escape from a meta-basin may be taking place 
already above T c . In this context, it is interesting to 
note that some theories of slow dynamics posit that ac- 
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FIG. 1: IS (solid line) and meta-basins (dotted line) visited 
during the course of a 70 particle MD run at T* — 0.49. The 
dotted line is slightly shifted along the energy axis for better 
clarity. IS transitions within a meta-basin are denoted /3 while 
transitions from one meta-basin to another are denoted a. 
The IS energies are divided by the total number of particles 
in the system. 



tivated processes are directly responsible for the onset of 
non-Arrhenius relaxation at temperatures far above T c 

s 

A simple way of rationalizing glassy dynamics in the 
regime where activation between attractors is the domi- 
nant mechanism is provided by the "trap" model advo- 
cated in [|l9|. In this picture, the dynamics of a coherent 
subregion of the system is summarized by the motion 
of a single point evolving in a landscape of "valleys" or 
"traps" , separated by barriers that can only be overcome 
via activation. The wandering of the phase point is de- 
scribed by a simple master equation with hopping rates 
that encode the statistics of barrier heights and the ge- 
ometry of phase space. The simplest of the trap models 
assumes that the top of the landscape is flat, and thus 
the rate of escape from a trap is related only to the trap 
depth. If one further assumes that the probability to 
reach any new trap is identical, the specification of the 
distribution of trap energies determines the behavior of 
all dynamical observables. Interestingly, Ben Arous et 
al. |2Cfl have shown that for a finite size p-spin model, 
the short time dynamics are identical to that predicted 
by mct, but the long time dynamics are precisely given 
by the trap model predictions. The purpose of this work 
is to make concrete connections between the landscape 
properties of a model glass forming liquid, and the trap 
picture above the glass temperature. The relevance of 
the trap picture to describe aging dynamics in the glass 
phase will be discussed in a later work. 

The system that we study is the modified 80-20 
Lennard- Jones system studied by Sastry et al. and 
we refer the reader to this work for the details. Since the 
macroscopic liquid should be viewed as a partitioning 



of noninteracting 'coherent' cells within which the land- 
scape picture is sensible ]|, |l^, |ll|, [2l|| , we study small 
systems containing between 70 and 150 particles. Long 
(up to 2 x 10 9 total time steps) series of inherent struc- 
tures are produced from long molecular dynamics runs at 
various temperatures between T* — 0.85 and T* = 0.49 
(T c w 0.45). This corresponds to a total time run of up to 
13 /is. Quenches are generated at intervals much shorter 
than the a-relaxation time at a given temperature from 
a combined steepest descent and conjugate gradient al- 
gorithm. This procedure is sufficient to resolve nearly all 
significant inherent structure transitions. 

In Fig.l we show a portion of an is history verses time 
(Eis(t) vs t) for T* = 0.490. One notes several remark- 
able features in this time series pT| . First, long regions 
where the system switches between only a few elemen- 
tary IS exist. Such long-lived collections of IS clearly 
provide deep trap states. Using the meta-basin defini- 
tion of Buchner and Heuer j^] , the time series of IS may 
be mapped into a time series of the visited meta-basins. 
Transitions between is within a meta-basin involve very 
small flexing of a cage, while transitions between meta- 
basins (signifying release from a trap) involve collective 
rearrangements of particles. This observation is fully con- 
sistent with Stillinger's landscape picture of the (3 and a 
processes, whereby the (3 process corresponds to tran- 
sitions between is within the meta-basin, while the a 
process corresponds to transitions between meta-basins 
Visual inspection demonstrates clearly that deeper 
trap states are longer lived, in accordance with the trap 
model [^9). It is important to note that although the is 
energy is used a label of the instantaneous state of the 
system, no assumptions are made concerning the transi- 
tion dynamics between meta-basins. In particular, non- 
trivial entropic factors may contribute to dynamics of 
meta-basin transitions. 

In simple trap models, the fundamental quantity is 
the distribution of trap depths. The most commonly 
used model assumes an exponential distribution of trap 
energies, i.e. p(E) = (1/T )exp(-E/T ). This model 
yields a power law distribution of trapping times ip(r) ~ 
—(1+—) 

t T o\ and a power law correlations for simple dynam- 
ical variables |l^, |2|, ^4| f25|. This model also has 
a strict glass transition at the temperature Tq. Other 
trap models, such as the Gaussian model with p(E) = 

( CXP ^ ^^wEo ^-^ j a l so display glassy phenomenology, 
such as a super- Arrhenius growth of the relaxation time, 
stretched exponential decay and (interrupted) aging jl9| . 

Assuming that the meta-basins form traps, the com- 
puted distribution of trap energies may be calculated di- 
rectly. At all temperatures studied, the distribution is 
well fitted by a Gaussian, as shown in Fig. 2(a). Previous 
studies have demonstrated that the distribution of inher- 
ent structures (not meta-basins) is Gaussian jl0|, ^6], p7| . 

Given simple activated dynamics to leave a trap of 
depth E, i.e. trap lifetimes r = Toexp(f3E), the dis- 
tribution of trapping times may be computed. For the 
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FIG. 2: (a) Distribution of meta-basin energies (per particle) 
for T* = 0.490, 0.575, 0.669 and 0.764 (from left to right). 
The solid lines are the Gaussian fits for the corresponding 
temperature, (b) Distribution of hopping times from meta- 
basins in a 70 particle system for T* = 0.764, 0.669, 0.575, 
and 0.490 (from left to right). The curves are shifted in time 
to fit inside the figure. 

Gaussian trap model, this yields a distribution of trap- 
ping times ip(r) = exp f- ^ ln W T ") ^ j ; wne re 

A = f3E . Note that the product Tip(r) is log-normal. 
In Fig. 2(b) we show fits of k^-n/^r)) vs. ln(r). The log- 
normal behavior predicted by the Gaussian trap model 
fits remarkably well over a wide range of time scales. 
The attempt time to is determined to be on the order 
of ~ lps, reasonable for the "vibrational" prefactor of an 
activated process in a dense Lennard- Jones liquid. Given 
that meta-basins can have lifetimes spanning up to a frac- 
tion of one /is, there clearly exists a separation of time 
scales between inter- and intra-trap (meta-basin) motion. 
This timescale separation does not exist at the level of 
single IS above T c . It is interesting to note that previous 
computational studies have determined power law wait- 
ing time distributions in a variety of systems |0, E8L E9| . 




FIG. 3: Correlation function C(t) as defined in the text. Both 
power law and stretched exponential fits are shown. The data 
clearly supports the stretched exponential fit over the power 
law fit. The solid line is for T* = 0.490 and dotted line is for 
T*=0.575. 



On the basis of a plot of ln-0(r) vs. m(r) it is diffi- 
cult to detect the distinction between the power law and 
log-normal versions of ip(r), while plots of 1b.(t0(t)) vs. 
ln(r) clearly reveal this distinction. The Gaussian widths 
A determined from fits of i/j(t) are larger than /3Eq de- 
termined from p(E) directly. This fact shows that one of 
the crude assumptions of the trap model cannot be taken 
literally; for example the energy depth of the meta-basin 
alone does not determine the absolute barrier for trap 
escape. Interestingly, the nontrivial entropic and saddle 
point dependence of Eq only quantitatively, and not qual- 
itatively affects estimates of the distribution of trapping 
times. 

We now compare simple dynamical predictions of the 
Gaussian trap model to computer simulations in the 
Lennard-Jones mixture. The basic quantity investigated 
is C(t) = (5eT(t)5eT(0)}, namely, the correlation func- 
tion of fluctuations of visited meta-basin depths. In the 
exponential trap model, C(t) decays as a power law, 
while for the Gaussian model the behavior of C(t) may 
be approximated with a stretched exponential decay |l9| : 

C(t) ~ exp (— (a(T)t b ( T })) , with a stretching coefficient 

_ i 

b(T) ~ (l + i(A) 2 ) 2 . For the range of temperatures 
studied, the Gaussian trap prediction of stretched ex- 
ponential relaxation is reasonably well borne out, as is 
shown in Fig. 3. Approximate stretched exponential de- 
pendence is seen for all temperatures T* < 1. At the two 
lowest temperatures studied in this work (T* = 0.575 
and T* = 0.490) where the extracted values of 6(T) are 
most accurate, the extracted values of b(T) are 6(0.575) 
= 0.68±0.1 and 6(0.490) = 0.61±0.05, respectively. The 
values predicted from the trap model are 6(0.575) = 
0.45±0.05 and 6(0.490) = 0.41±0.1, respectively. Thus, 
the predicted values are in reasonable agreement with the 
values extracted from simulation |30[ . 

Our results show that activated process are important 
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even above T c in supercooled liquids. This does does not 
contradict the applicability of mct to describe short time 
dynamics in this regime, nor does it contradict the results 



of ItS, 
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which show that most accessible saddles 
are unstable above T c . The long time dynamics, how- 
ever, will be dominated by rare, deep traps. Our work 
is consistent with this, provided that traps (meta-basins) 
smoothly transition from a collection of connected inher- 
ent structures above T c to single inherent structures be- 
low T c . That is, above T c the barriers separating the in- 
herent structures inside a meta-basin are small compared 



to temperature, while below T c they are large. Com- 
puter simulation evidence supports this picture fl6|| . This 
would imply that the simple "single-level" trap model 
above T c should become a "multi-level" trap picture be- 
low T c ||, where interestin g a ging effects should then 
take place as in spin-glasses fl32f . 
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